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Abstract 


A  finite  element  investigation  was  conducted  to  analyze 
a  crack  growing  from  a  round  (large  radius)  notch  under 
cyclic  conditions  at  1200 °F.  The  results  of  the  finite  ele¬ 
ment  analysis  were  compared  to  previously  published  results 
obtained  for  a  crack  emanating  from  a  blunt  notch  and  a 
single-edge-cracked  specimen. 

The  finite  element  program  used  was  a  two-dimensional 
materially  and  geometrically  nonlinear  finite  element  code 
called  SNAP.  The  program  has  the  capability  to  release 
fixed  nodes  to  allow  crack  growth  and  the  ability  to  simu¬ 
late  crack  closure  under  cyclic  conditions.  Constitutive 
equations  set  forth  by  Bodner  and  Partom  were  used  to 
account  for  the  nonlinear ,  viscoplastic  material  behavior 
exhibited  by  IN-718  at  1200 °F,  The  load  spectrum  included 
loading  under  a  positive  load  range  (R  ■  0.1)  and  loading 
under  a  negative  load  range  (R  *  -1.0). 

Finite  element  analysis  of  a  crack  growing  from  a  round 
notch  under  cyclic  loading  provided  cr*ck  opening  profile 
information*  opening*  and  closing  loads*  stress  and  strain 
profiles  and  plastic  zone  estimations. 

Crack  closure  develops  over  a  longer  distance  for  the 
round  notch  than  that  for  a  blunt  notch.  The  elastic  stress 
intensity  influence  on  the  round  notch  was  found  to  be 
within  one  radius  of  the  notch.  In  general*  the  crack 
growth  trends  observed  for  the  round  notch  agree  with  those 
seen  in  the  blunt  notch  and  single-edge-crack  specimen. 
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I.  Introduction 


Background 

Aircraft  engine  components  and  structures  often  operate 
under  conditions  severe  enough  that  the  designer  must  ensure 
that  the  component  or  structure  has  an  adequate  fatigue 
life.  The  United  States  Air  Force  has  taken  steps  to  ensure 
that  fatigue  resistance  is  considered  a  primary  design  cri¬ 
terion  by  developing  the  Engine  Structural  Integrity  Program 
(ENSIP).  This  program  provides  an  organized  approach  to  the 
design,  analysis,  and  verification  testing  of  all  newly 
developed  aircraft  engines  (1). 

The  ENSIP  assumes  that  flaws  exist  and  may  be  found  by 
means  of  inspection  or  by  crack  growth  analysis  if  the  flaw 
is  undetectable  by  nondestructive  inspection  (NDI)  tech¬ 
niques.  This  is  by  no  means  straightforward,  since  crack 
growth  analysis  comprises  both  crack  initiation  as  well  as 
crack  propagation.  For  simple  geometries,  linear  elastic 
fracture  mechanics  (LEFM)  provide  a  good  estimate  of 
remaining  useful  component  life.  However,  for  components  in 
the  latter  stages  of  the  compressor  and  turbine,  tem¬ 
peratures  may  reach  1400 °F.  Thus,  crack  propagation  will  be 
influenced  by  the  superposition  of  high  temperatures  and 
stresses  and  LEFM  crack  growth  analysis  will  not  provide 
accurate  engine  component  life  determination. 

Most  fatigue  failures  originate  at  some  kind  of  stress 
concentration  or  notch.  At  the  roots  of  notches  the  stress 
may  exceed  the  yield  stress  of  the  material.  The  growth  of 
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small  cracks  (small  compared  to  the  notch  dimensions)  in 
such  regions  will  be  controlled  by  local  notch  plasticity. 
When  the  crack  is  propagating  in  the  notch  plastic  field  it 
becomes  incorrect  to  use  LEFM  to  characterize  crack  growth. 
Thus,  elastic-plastic  or  viscoplastic  analysis  techniques 
must  be  used.  This  study  used  the  finite  element  method  to 
evaluate  crack  growth  at  a  round  notch  using  a  viscoplastic 
constitutive  law. 

Approach 

As  previously  stated,  this  research  concentrated  on  a 
crack  emanating  from  a  round  notched  compact  tension  speci¬ 
men.  The  specimen  studied  was  IN-718,  a  nickel-based 
superalloy  used  in  the  F-1QG  engine.  SNAP,  a  nonlinear 
finite  element  program  developed  by  Brockman  (2)  was  used 
for  this  analysis.  This  finite  element  program  is  a  two- 
dimensional  plane  stress  and  plane  strain  code  which 
includes  material  and  geometric  nonlinearities.  For  this 
study,  eight-noded  quadratic  isoparametric  elements  were 
used,  SNAP  has  the  ability  to  release  fixed  nodes  to  simu¬ 
late  crack  growth  and  also  simulate  crack  closure  under 
cyclic  conditions.  In  addition,  a  set  of  constitutive 
equations  called  the  Bodner-Partom  viscoplastic  flow  law  (3) 
was  used  to  model  the  material  behavior  during  load  incre¬ 
ments.  Plastic  strains  were  calculated  using  an  implicit 
technique  developed  by  Brockman  (2)  and  were  incorporated 
into  the  finite  element  model  by  the  residual  force  method. 
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The  primary  goal  of  the  study  was  to  analyze  a  crack 
growing  from  a  round  notch  specimen  and  compare  these 
results  to  those  obtained  for  a  blunt  notch  (4)  and  a 

s 

single -edge -cracked  specimen  (5).  For  comparison  to  the 
blunt  notch  specimen,  load  was  input  as  a  saw-toothed 
stress-time  pattern  of  constant  amplitude  with  a  load  ratio 
of  0.1.  Maximum  load  amplitudes  were  applied  to  result  in 
two  load  levels  for  which  the  notch  was  elastic  and  plastic. 
Cyclic  load  frequency  was  maintained  at  1.0  Hz.  For  com¬ 
parison  to  the  single-edge-cracked  specimen,  the  loading  was 
of  the  same  pattern  but  with  a  load  ratio  of  -1.0  and  for 
the  plastic  loading  case  only. 

Literature  Review 

As  stated  earlier,  the  problem  of  short  crack  growth  at 
notches  is  important  because  most  fatigue  failures  originate 
at  some  kind  of  stress  concentration  and  the  period  of  crack 
growth  through  the  notch  stress  field  may  represent  a  major 
part  of  the  total  life.  Thus,  a  proper  knowledge  of  short 
crack  growth  from  notches  is  essential  for  accurate  predic¬ 
tions  of  the  fatigue  life. 

Short  cracks  tend  to  grow  faster  than  predicted  using 
LEFM  principles.  Leis  (6)  observed  that  LEFM  is  limited  to 
small  scale  yielding  at  the  crack  tip.  The  size  of  the 
plastic  zone  is  usually  much  smaller  than  the  crack  length. 
The  author  found  that  short  cracks  have  large  plastic  zones 
that  are  generally  the  same  size  as  the  crack  length.  As  a 
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consequence,  the  use  of  LEFM  leads  to  unconservative  esti¬ 
mates  of  the  short  crack  fatigue  life.  Smith  and  Miller  (7) 
observed  that  when  a  crack  emanates  from  the  notch  root,  it 
first  grows  under  notch  plasticity  control  and  at  a 
decreasing  rate  because  of  the  falling  notch  strain  field. 

As  the  crack  reaches  the  end  of  the  notch  plastic  zone,  it 
then  grows  under  its  own  plasticity  control. 

Several  attempts  have  been  made  to  account  for  the 
faster  growth  of  cracks  at  notches.  In  general,  most  pro¬ 
vide  corrections  to  the  physical  crack  length  to  account  for 
the  effect  of  the  notch.  El  Haddad  et  al.  (8)  modified  the 
stress  intensity  range  (AK)  by  adding  an  artificial  length, 
1Q,  to  the  crack  length,  a: 

Ak  =  Ao  /Or  (a  +  10)  )l/2  (1.1) 

where  Acr  is  the  local  stress  range.  For  elastic-plastic 
solutions,  El  Haddad  et  al.  (9)  proposed  using  either  the 
J-integral  range  or  a  strain-based  intensity  factor,  c, 
defined  as. 

AKc  *  EAe/Uta  +  l0))l/2  (1.2) 

where  E  is  the  elastic  modulus,  a e  is  the  local  plastic 
strain  range,  and  (a  +  1Q)  is  the  corrected  crack  length. 
This  formulation  gave  acceptable  results  for  notched  and 
unnotched  geometries. 

In  recent  years,  fatigue  crack  closure,  originated  by 
Elber  (10),  has  been  used  to  explain  the  behavior  of  cracks 
growing  at  notches.  Elber  found  that  crack  closure  occurred 


15 


as  a  consequence  of  crack  tip  plasticity.  At  the  tip  of  a 
growing  fatigue  crack,  a  zone  of  residual  deformation, 
caused  by  the  plastic  zone  during  increased  loading,  is  left 
in  the  wake  of  the  advancing  crack  tip.  When  the  fatigue 
load  decreases,  this  residual  deformation  causes  the  crack 
to  close.  According  to  this  concept,  the  crack  cannot  pro¬ 
pagate  while  it  is  closed.  Thus  AK  is  replaced  by  AKeff 
and  is  defined  as: 

AKeff  =  Kmax  “  Kcl  (1.3) 

where  K^x  is  the  maximum  stress  intensity  and  Kc^  is  the 
stress  intensity  value  for  which  the  crack  tip  closes  (or 
opens)  during  the  loading  cycle.  Newman  (11)  performed  an 
elastic-plastic  finite  element  analysis  of  cracks  growing 
from  notches.  He  used  spring  elements  along  the  crack  to 
satisfy  changing  boundary  conditions.  He  found  that  Kc^  can 
explain,  in  terms  of  AKeff,  the  non-LEFM  behavior  of  short 
cracks  growing  at  notches.  Nicholas,  Palazotto  and  Bednarz 
(5)  performed  a  numerical  analysis  which  modeled  the  plasti¬ 
city  induced  closure  of  a  short  crack  emanating  from  a 
single-edge-cracked  specimen.  They  used  a  two-dimensional 
plane  stress/plane  strain  viscoplastic  finite  element  code 
callod  VISCO,  which  was  originally  developed  by  Hinnerichs 
(12).  They  examined  crack  growth  and  crack  closure  for 
TI-6246  at  room  temperature  for  stress  values  of  approxima¬ 
tely  .60  and  .90  material  yield  strength  at  load  ratios  of 
-1.0  and  0.1.  Finite  element  results  included  displacements 
along  the  crack  surface,  stress-strain  values  immediately  in 
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front  of  the  crack  tip,  stress  profiles  ahead  of  the  crack 
tip,  plastic  strains  behind,  as  well  as  in  front  of,  the 
crack  tip  and  crack  closure  information. 

Mercer  (4)  performed  a  viscoplastic  finite  element  anal¬ 
ysis  of  a  crack  growing  from  a  blunt  notch  under  cyclic 
loading  conditions  at  1200 °F.  The  two-dimensional  finite 
element  code  called  SNAP  was  modified  by  Mercer  to  include 
the  capability  to  simulate  crack  closure  by  the  addition  of 
springs  to  the  nodes  along  the  crack  line.  This  method  is 
similar  to  the  one  used  by  Newman  (11).  He  examined  crack 
growth  and  crack  closure  for  IN-718  at  elastic  and  plastic 
loading  for  a  load  ratio  of  0.1.  The  finite  element  analy¬ 
sis  provided  crack  tip  stress  and  strain,  plasticity  induced 
closure  effects  and  crack  opening  profile  information. 

The  review  of  the  literature  indicates  that  the  behavior 
of  cracks  is  not  completely  understood  and  additional 
research  needs  to  be  done.  In  addition,  numerical  methods, 
especially  finite  element  analyses,  can  be  used  to  provide 
helpful  information  in  predicting  the  behavior  of  cracks  at 
notches. 
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II.  Theory 


Linear  Elastic  Fracture  Mechanics  (LEPM) 

Linear  elastic  fracture  mechanics  relates  the  stress 
field  and  displacements  near  a  crack  tip  to  the  stress 
applied  to  the  component  for  various  crack  geometries.  The 
difference  between  one  cracked  component  and  another  lies  in 
the  magnitude  of  the  stress  field  parameter,  Kj,  defined  as 
the  mode  I  (opening  or  tensile  mode)  stress  intensity  fac¬ 
tor.  In  essence,  the  whole  stress  field  at  the  crack  tip  is 
known  when  the  stress  intensity  solution  i3  known  (13).  The 
stress  intensity  parameter  has  been  discussed  explicitly  in 
fracture  mechanics  literature  (13,14),  and  thus  will  not  be 
elaborated  upon  herein.  However,  the  following  stress 
intensity  factors  for  some  of  the  more  well-known  specimen 
geometries  are  given: 

Center  Cracked  Plate  (14) : 

Kj  *  o(«a  )  1/2 (aec  Tr_a)l/2  (2.1) 

W 

where  o  is  the  notch  stress ,  a  is  the  crack  length  measured 
from  edge  of  notch,  and  W  is  the  specimen  width. 

Single  Edge  Notched  Plate  (13): 

K i  »  1.12  c(* a)!/2  (2.2) 

where  o  is  the  notch  stress  and  a  is  the  crack  length 
measured  from  edge  of  notch.  This  equation  is  only  valid 
for  small  cracks. 
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Compact  Tension  Specimen  (15) : 


Kj  =  _P _ .  f  (a)  (2.3) 

BWl/2  w 

where 

f  (a)=  (2  +  a  )  [0.866  +  4.64a  -  13.32a2  +  14.72  a3  -  5.6a4] 

W  (1  -a) 3/2 

P  is  the  applied  load,  B  the  specimen  thickness,  W  the  spe¬ 
cimen  width  and  a  =  A  where  A  is  measured  from  point  of 

W 

load  application.  This  solution  is  valid  for  A  >  0.2. 

w  “ 

Bodner-Part om  Constitutive  Law 

The  Bodner-Partom  flow  law  accounts  for  viscoplastic 
behavior  as  well  as  rate  sensitivity  and  strain  hardening 
effects  (3).  Work  has  been  previously  presented  on  the 
ramifications  of  this  law  (16),  but  for  completeness,  basic 
concepts  and  equations  are  presented.  After  initial 
yielding  the  material  behavior  will  be  partly  elastic  and 
partly  plastic.  During  any  increment  of  stress,  the  changes 
of  strain  are  assumed  to  be  separable  into  elastic  and 
plastic  components  such  that 

*  «  6  •  p 

eij“cij+eij  (2.4) 

The  plastic  strain  rate  can  be  expressed  in  the  form  of 
the  Prandtl-Reuss  equation  of  classical  plasticity  as 

cp  x 

E i j  ”  X  sij  (2-5) 
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where  X  is  the  proportionality  constant  called  plastic 
multiplier  and  S^j  is  the  deviatoric  stress  tensor.  The 
form  of  X  depends  on  the  flow  law  used.  Bodner  expressed 
X  as  a  function  of  the  deviatoric  stress  tensor  and  the 
plastic  strain  rate  tensor  and  can  be  written  as 


X2 


P 

»  p2 

J2 


(2.6) 


P  *P  -P 

where  D2  =  1/2  cfjeij  and  is  called  the  second  invariant 
of  the  plastic  strain  rate  tensor  and  J2  *  1/2  S^jS^j  is 
called  the  second  invariant  of  the  deviatoric  stress  tensor. 
Bodner  and  Partom  expressed  D2P  as 

n 


P  2 

D2  =  %  exp 


re)  m 


(2.7) 


where  0o  is  the  limiting  value  of  plastic  strain  rate  in 
shear,  Z  is  the  measure  of  material  hardening,  and  n  is  a 
rate  sensitivity  parameter.  Z  depends  on  the  deformation 
history  of  the  material  and  is  assumed  to  be  a  function  of 
the  plastic  work,  Wp,  such  that 

Z  -  Zi  +  (Zq  -  Zi)  exp  (-  mwA  (2.8) 

\  w 

where  Zi  and  20  are  the  material's  maximum  and  initial  value 
of  hardness  respectively,  and  m  is  a  constant  that  controls 
the  rate  of  work  hardening. 


Generally,  the  plastic  work  is  expressed  as 

•  P 

Wp  =| Si j  dt 


(2.9) 
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However,  to  account  for  the  thermal  recovery  of  hardening  at 
high  temperature,  the  plastic  work  must  be  redefined  as 


Wp  = 


lsij  Cijdt  +  I  Zrec  dt 
m(Zi  -  Z) 


(2.10) 


where  the  rate  of  thermal  recovery  is  defined  as 

) r 

Zrec  =  -  AZ]l  Zi  ~  J  (2.11) 

and  Z2  is  the  value  of  Z  that  gives  the  minimum  expected 
value  of  hardening  at  a  given  temperature.  A  and  r  are 
material  constants  chosen  to  match  creep  test  data. 

The  constants  used  for  IN  718  as  determined  by  Beaman 
(17)  are: 


Dq  =  106  sec-1 
Z0  =  235.3  Ksi 
Zi  *  260.3  Ksi 
Z2  *  104.1  Ksi 
n  “3.0 

m  «  2.875  Ksi”1 
A  “  1.5  X  10”1  sec-1 


r  *  7.0 

Viscoplastic  Solut ion  Procedure 
Bodner-Partom  Solution 

SNAP  employs  a  semi -implicit  technique  to  solve  the 
Bodner  equations.  Plastic  strains  are  calculated  based  upon 
a  semi-implicit  scheme  which  averages  the  strain  rates  at 
the  beginning  and  end  of  each  increment  (2).  Initial  esti¬ 
mates  of  conditions  at  the  end  of  the  increment  are  obtained 
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using  Euler  extrapolation.  This  scheme  is  summarized  as 
follows : 

1.  Initial  values  for  stress,  plastic  strain,  plastic  work, 
Z,  and  the  strain  increment  for  the  load  step  are  input  into 
the  Bodner  routine. 

2.  The  following  Bodner -Par tom  equations  are  solved  for 
this  increment: 


(2.12) 


(2.13) 

(2.14) 


(D2P)iy/2 


^2 


i-1 


(2.15) 


(2.16) 


(2.17) 


i-1 


Zrec^ 


Wp 


1  «  Wpi-i  +  (SijC  +  m(Zi  -  Z1)  (2. 

3.  New  estimates  of  stresses,  plastic  strains,  plastic 
work,  and  Z  are  calculated  using  an  average  of  rates.  As 


18) 
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described  by  Owen  and  Hinton  (18),  the  plastic  strain  rate 


is: 


i 

=  dti 


(1**0) 


i-1 


where  in  the  semi-implicit  case  a  =  1/2. 


(2.19) 


4.  The  difference  in  stress  values  from  steps  (2)  and  (3) 
are  calculated.  If  the  difference  is  within  a  prescribed 
tolerance,  the  solution  for  this  increment  is  complete  and 
the  next  increment  begins.  If’  the  difference  is  not  within 
the  allowable  tolerance,  a  convergence  loop  is  entered. 


Solu ti on  P rocedure 

The  overall  solution  procedure  employed  in  SNAP  is  the 
residual  force  method.  The  residual  force  method  uses  the 
elastic  stiffness  matrix  during  the  entire  analysis  and 
incorporates  plasticity  effects  by  including  a  plastic  load 
vector.  The  residual  force  method  has  the  following  form: 

*  [fji  +  {i}1"1  (2.20) 

where  (K]  is  the  constant  elastic  stiffness  matrix,  ^  are 
the  nodal  displacements  vector,  |t^ *  is  the  applied  load 
vector,  and  ^ is  the  element  internal  force  vector 
developed  by  the  accumulation  of  plastic  deformation.  The 
superscript  i  reflects  the  current  increment.  The  visco¬ 
plastic  solution  for  a  typical  time  step  using  the  residual 
force  method  proceeds  as  follows: 

1.  Compute  the  current  time 

ti  «  dti  +  ti’1  (2.21) 
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2.  Calculate  the  plastic  strain  rate  using  the  Bodner  model 
relations 


(2.22) 


3.  Compute  the  plastic  strain  using  the  semi-implicit 
scheme 


where  e  =  1/2  for  the  semi -implicit  technique.  For  0  »  0, 
we  get  the  Euler  integration  scheme  (fully  explicit)  and  for 
G=  1  we  get  the  fully  implicit  scheme. 

4.  The  plastic  load  vector  is  computed  by 


[BHD] 


VOL 


i 

dVOL 


(2.24) 


where  [Bl  is  the  strain-displacement  matrix  and  tDl  is  the 
elastic  material  property  matrix. 

5.  Compute  the  current  external  load  vector 


Qi  at*  ♦  {t}1-1 


(2.25) 


where  T  is  the  known  force  rate  vector. 

6.  Compute  the  nodal  displacements  using  Gaussian  elimina¬ 
tion  techniques 

{u}  1  *  IK  pi  (  £r}i  +  {i)1’1)  (2.20 

7.  Compute  the  total  strain  from  the  strain-displacement 
relationship 
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8. 


Compute  the  current  stress 


9.  The  particular  time  step  is  determined  when  the  differ¬ 
ences  in  stress  values  discussed  earlier  are  within  a  set 
tolerance.  If  convergence  is  not  obtained  within  40  itera¬ 
tions,  the  time  step  is  halved  and  steps  1-9  are  repeated 
for  this  particular  time  step.  If,  after  10  successive  time 
step  reductions,  convergence  is  not  achieved,  the  problem 
solution  is  stopped. 


Crack  Closure  and  Gr owth 

Crack  growth  and  closure  were  implemented  by  Mercer  (4) 
into  the  SNAP  program  by  the  addition  of  "springs*  to  the 
boundary  node  along  the  crack  line.  This  method  is  similar 
to  that  used  by  Newman  (19).  The  spring  stiffness  for  each 
node  was  related  to  the  boundary  condition  of  that  respec¬ 
tive  node.  If  the  spring  was  made  very  stiff  the  node  was 
considered  fixed  and  if  the  node  was  free  the  stiffness  was 
set  equal  to  zero. 

Por  the  nodes  lying  along  the  crack  line,  a  flag  for 
each  node  behind  and  ahead  of  the  crack  tip  was  monitored. 

A  positive  value  indicated  a  fixed  node  ar.d  a  negative  value 
meant  the  node  was  free.  During  each  load  increment,  the 
nodal  displacements  along  the  crack  line  were  also  monitored 
to  determine  whether  the  node  had  closed  (negative  displace¬ 
ment)  or  opened  (positive  displacement).  If  the  node  had 
closed,  the  flag  was  set  equal  to  a  positive  value  and  the 
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spring  stiffness  was  set  to  1000  times  the  maximum  diagonal 
entry  in  the  structure  stiffness  matrix  and  the  stiffness 
matrix  was  updated  to  account  for  the  changing  boundary  con¬ 
ditions.  If  the  node  had  opened,  the  flag  was  set  equal  to 
a  negative  value  and  the  spring  stiffness  was  set  equal  to 
zero  and  the  stiffness  matrix  was  recalculated.  This  proce¬ 
dure  of  adding  or  deleting  spring  stiffnesses  to  the  struc¬ 
tural  stiffness  matrix  had  the  advantage  that  the  global 
stiffness  matrix  did  not  have  to  be  reformulated  and  decom¬ 
posed  whenever  the  crack  opened  or  closed  (4). 

Crack  growth  was  accomplished  by  releasing  nodes  when  a 
user-defined  criterion  such  as  stress  or  load  was  exceeded. 
The  selected  method  released  an  element  with  every  cycle  of 
loading  at  98  percent  load  for  R  =  0.1,  where  R  is  the  mini¬ 
mum  to  maximum  load  ratio.  For  R  =  -1.0,  an  element  was 
released  with  every  cycle  of  loading  at  98  percent  of  the 
stress.  These  methods  allowed  for  full  recovery  of  crack 
tip  plastic  strain  in  each  cycle,  and  no  convergence 
problems  were  encountered  (4). 

J-Integral 

The  elastic  stress  intensity  factor  can  be  determined 
for  crack  problems  by  using  finite  element  methods  including 
the  J-integral,  compliances,  and  extrapolation  of  the  stress 
and  displacement  fields  near  the  crack.  The  J-integral 
technique  seems  to  provide  the  best  results  (4);  thus,  the 
J-integral  concept  was  used  in  this  study. 
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The  J-integral  concept  is  baaed  upon  an  energy  balance 
approach.  The  concept  was  introduced  by  Rice  (20).  For  an 
elastic  case,  J  is  equivalent  to  the  energy  release  rate  G 
(13).  The  J-integral  formula  represents  an  energy  influx 
across  an  arbitrary  boundary  r  around  the  crack  tip  and  is 
defined  as: 

J  =  [Wdy  -  Ti  3Ui  ds] 

r  3  x  (2.29) 

where  W  is  the  strain  energy  density,  T^  is  the  traction 
vector,  and  is  the  displacement  vector  along  the  contour. 
Obtaining  solutions  for  the  J-integral  in  actual  specimens 
is  difficult.  It  is  generally  necessary  to  use  finite  ele¬ 
ment  techniques.  Mercer  (4)  describes  in  detail  how  the 
actual  J-integral  numerics  were  computed  using  the  SNAP 
finite  element  program. 
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III.  Results  and  Discussion 


A  viscoplastic  analysis  of  a  crack  growing  from  a  com¬ 
pact  tension  specimen  with  a  round  notch  was  investigated 
using  the  SNAP  finite  element  program.  A  specimen  con¬ 
sisting  of  IN-718  at  1200°  was  modeled  under  cyclic  con¬ 
ditions  at  two  load  levels  with  a  load  ratio  of  0.1.  The 
specimen  was  also  analyzed  under  plastic  loading  with  a  load 
ratio  of  -1.0.  Both  elastic  and  viscoplastic  analyses  were 
performed  on  an  uncracked  specimen  considering  a  plane 
stress  approach.  A  viscoplastic  analysis  was  performed  for 
a  crack  growing  from  an  initial  length  of  0.0  inches  to  a 
final  crack  length  of  0.024  inches. 

Data  obtained  from  the  finite  element  model  analysis 
included  an  elastic  stress  intensity  solution,  crack  opening 
displacements,  opening  and  closure  loads,  stress  values  in 
front  of  the  crack  tip,  and  plastic  strains  along,  as  well 
as  in  front  of,  the  crack  tip.  These  values  are  compared  to 
data  obtained  for  a  crack  emanating  from  a  blunt  notch  and  a 
single -edge-cracked  specimen. 

Compact  Tension  Finite  Element  Modeling 

The  compact  tension  specimen  had  a  round  notch  with  a 
radius  of  .197  inches.  The  specimen  geometry  is  shown  in 
Figure  3.1.  Due  to  symmetry,  only  half  of  the  specimen  was 
modeled.  The  model  used  the  eight-noded  quadratic  isopara¬ 
metric  elements. 
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The  overall  finite  element  model  is  shown  in  Figure  3.2. 
The  refined  elements  used  in  the  neighborhood  of  the  crack, 
along  with  a  typical  Gauss  point  location,  are  shown  in 
Figure  3.3.  This  refinement  is  identical  to  that  used  by 
Mercer  (4)  for  a  blunt  notch  (r  =  0.47  in).  The  number  of 
elements  near  the  crack  tip  and  the  degree  of  model  refine¬ 
ment  represented  a  trade-off  between  the  amount  of  crack 
growth  required  and  the  computer  storage  space  available,  as 
well  as  the  computer  time  required  to  carry  out  the  simula¬ 
tion.  The  model  with  crack  tip  refinements  had  about 
1000  degrees  of  freedom,  490  nodes,  and  140  elements.  The 
crack  tip  elements  were  0.002  inch  square. 

Mercer  (4)  performed  an  extensive  series  of  validation 
tests  using  the  SNAP  program.  These  tests  include: 
a)  V-notched  elastic-plastic  tension  problem  as  described  by 
Yamada  (21),  b)  center-cracked  panel  analysis  as  reported  by 
Hinnerichs  (12)  and  Henkel  (22),  and  c)  compact  tension 
specimen  analysis  as  examined  by  Wilson  (23).  The  SNAP 
program  was  found  to  provide  satisfactory  results  for  all  of 
the  above  problems.  Since  acceptable  accuracy  was  shown  by 
Mercer  using  the  SNAP  program,  no  extensive  accuracy  tests 
were  carried  out  for  the  round  notch  study.  However,  a 
validation  run  was  performed  for  the  round  notch  uncracked 
specimen  and  SNAP  gave  acceptable  results  as  will  be  seen  in 
the  following  section. 
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Uncracked  Specimen  Elastic  Analysis 


Irregularities  in  a  component/  such  as  holes  and 
notches,  will  have  high  localized  stresses  and  may  also  pro¬ 
vide  possible  sites  for  fatigue  in  a  cyclic-loaded  com¬ 
ponent.  These  high  localized  stresses  are  called  stress 
concentrations  and  they  provide  useful  information  regarding 
the  irregularity’s  stress  distribution  and  are  usually 
expressed  in  terms  of  a  stress  concentration  factor.  This 
factor  was  determined  for  the  round  notch  condition  using 
the  finite  element  method  and  was  then  compared  to  pre¬ 
viously  published  compact  tension  stress  concentration 
results. 

The  finite  element  model  was  monotonically  loaded  to  0.5 
Kips  using  linear  elastic  material  properties.  The  elastic 
stress  concentration  factor,  k,  is  defined  as  the  ratio 
between  the  maximum  stress  and  the  nominal  stress  in  a 
notched  sample  and  can  be  calculated  by: 

k  *  Oy  (3.1) 

on 

where  ^  is  the  calculated  stress  at  the  notch  and  o  n  is  the 
applied  nominal  stress  at  the  notch.  The  nominal  stress  for 
the  specimen  considered  in  Figure  3.1  is  defined  by  Wilson 
(24)  as: 

°n-P  1  pl+3(W+  An)“ 

B  ( W-An)  (_  Tw  -  An)_ 

where  P  is  the  applied  load,  B  is  the  specimen  thickness.  W 
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is  the  specimen  width/  and  An  is  the  distance  from  load 
point  to  notch  location.  The  calculated  stress  concentra¬ 
tion  factor  for  the  round  notch  was  2.18.  As  shown  in 
Figure  3.4,  the  round  notch  stress  concentration  factor 
result  agrees  with  that  reported  by  Wilson  (24).  Wilson's 
results  are  for  a  compact  tension  specimen  with  three  dif¬ 
ferent  notch  shapes.  Figure  3.4  is  drawn  to  show  how  the 
round  notch  result  compares  to  Wilson's  notch  shapes.  The 
symbol/  R,  is  a  nondimens ional  notch  radius  and  is  defined 
as: 

R  =  r _ 

W-An  (3.3) 

where  r  is  the  notch  radius,  W  is  the  specimen  width  and  An 
is  the  distance  from  load  point  to  notch  location.  The 
stress  concentration  factor  obtained  by  Mercer  (4)  for  a 
blunt  notch  (r  *  .047")  was  4.29.  The  difference  between 
the  k  values  is  obvious;  the  round  notch  has  a  lower  maximum 
notch  stress  than  that  of  a  blunt  notch. 

Uncracked  Specimen  Viscoplastic  Ana lysis 

In  order  to  determine  which  loads  will  be  used  in  this 
study,  we  must  first  determine  if  the  load  to  be  applied  is 
an  elastic  or  plastic  load.  This  can  be  accomplished  by 
evaluating  the  stress  concentration  factor  and  by  observing 
if  the  applied  load  produced  any  plastic  strain.  Both  of 
these  methods  can  be  evaluated  using  the  finite  element 
model. 
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The  finite  element  model  was  monotonically  loaded  to  2.0 
Kips  in  one  second  using  the  previously  described  Bodner 
model.  The  stress  concentration  factor  is  defined  by 
equation  (3.1).  A  strain  concentration  factor,  as  described 
by  Wilson  (24),  can  be  defined  as: 


where  E  is  the  elastic  modulus,  Ey  is  the  calculated  strain 
obtained  from  the  finite  element  program,  and  °n  is  the 
nominal  stress  from  equation  (3.2).  If  the  stresses  and 
strains  are  elastic,  the  stress  concentration  factor  and  the 
strain  concentration  factor  will  be  equal.  The  stress  and 
strain  concentration  factors  are  shown  in  Figure  3.5  as  a 
function  of  the  load  level.  The  stress  and  strain  con¬ 
centration  are  the  same  up  to  loads  equal  to  1.2  Kips,  thus 
there  was  no  plastic  deformation  at  the  notch  up  to  this 
load  level.  However,  for  loads  greater  than  1.2  Kips,  the 
effect  of  plasticity  at  the  notch  is  observed.  The  strain 
intensity  increases  due  to  the  plastic  strain  increasing, 
while  the  stress  intensity  decreases  since  the  notch  is  no 
longer  elastic. 

The  plastic  strain  at  each  element  Gauss  point  was 
observed  and  the  0.1%  plastic  strain  location  obtained  by 
interpolation.  This  plastic  strain  value  was  the  same  value 

used  by  Mercer  (4).  The  size  of  the  plastic  zone  as  defined 
by  eyP-  0.1%  is  shown  in  Figure  3.6.  The  loads  selected  for 
crack  growth  analysis,  at  a  load  ratio  of  0.1,  were  1.148 
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Sips  and  1.55  Kips.  We  see  that  for  the  lower  load,  there 
is  no  notch  plasticity  while  at  the  higher  load  the  notch 
plasticity  was  0.013  inches.  The  load  selection  criteria 
used  for  this  study  were  similar  to  those  used  by  Mercer 
(4).  The  lower  load  was  an  elastic  load,  while  the  higher 
load  produced  significant  notch  plasticity.  It  should  be 
noted  that  the  lower  load  chosen  for  tnis  study  coincided 
with  Mercer’s  plastic  load  for  the  blunt  notch  study.  By 
choosing  these  two  load  levels,  crack  growth  could  be 
observed  under  two  different  cases:  one  within  the  elastic 
notch  influence  and  one  within  the  plastic  notch  influence. 

The  model  was  subjected  to  a  cyclic  load  condition  to 
examine  the  response  at  the  notch  for  these  selected  load 
levels.  The  load  ratio  used  was  R  =  0.0.  At  P  *  1.148 
Kips,  we  s?e  that  the  material  stays  within  the  elastic 
regime  as  shown  in  Figure  3.7.  However,  at  P  =  1.55  Kips, 
we  see  plastic  deformation  during  the  initial  cycle,  but 
subsequent  cycle  response  was  elastic.  The  stress  profiles 
ahead  of  the  notch  at  100%  and  vtero  load  are  shown  in  Figure 
3.8  and  Figure  3.9  for  the  elastic  and  plastic  load  levels, 
respectively.  We  can  see  that  for  maximum  loading,  the 
material  starts  out  in  tension  near  the  notch  and  goes  into 
compression  away  from  the  notch.  This  observation  is  due  to 
the  "bending"  effect  taking  place  as  the  compact  tension 
specimen  is  being  loaded  as  shown  diagrammatically  in  the 
inset  in  Figures  3.8  and  3.9.  For  the  plastic  load,  we  see 
that  the  material  goes  into  compression  when  the  load  is 


33 


reduced.  This  observation  tells  us  that  a  local  negative  R 
ratio  is  seen  for  an  applied  load  ratio  of  0.0. 

Elastic  Stress  Intensity  Calculation 

As  discussed  earlier,  the  elastic  stress  intensity  is  a 
commonly  used  parameter  to  characterize  linear  elastic  frac¬ 
ture  mechanics.  The  stress  intensity  factor  incorporates 
geometrical  terms  as  well  as  the  stress  level.  The  elastic 
stress  intensity  for  a  crack  growing  from  a  round  notch  was 
determined  from  the  finite  element  model.  The  finite  ele¬ 
ment  results  were  compared  with  the  compact  tension  specimen 
long  crack  solution  to  determine  the  effect  of  notch 
influence.  The  solutions  are  based  on  an  elastic  load  of 
1.148  Kips  and  the  elastic  constitutive  model  in  SNAP  was 
used. 

First,  we  must  establish  bounds  on  the  elastic  stress 
intensity  solution.  This  was  established  by  using  the  com¬ 
pact  tension  solution  for  long  cracks  and  edge  crack 
approximations  for  short  cracks.  The  compact  tension  solu¬ 
tion  is  given  by  earlier  equation  (2.3)  from  ASTM  standard 

E-647-83  (15).  This  equation  is  valid  for  b  >  0.2.  This 

W  ~ 

solution  is  shown  as  the  slightly  inclined  curve  in  Figure 
3.10.  A  small  crack  in  a  notch  acts  like  an  edge  crack  in  a 
plate  under  a  far  field  stress  equal  to  the  concentrated 
stress  at  the  notch.  The  edge  crack  solution  for  small 
cracks  is  defined  by  equation  (2.2).  Observing  that  the 
magnitude  of  the  stress  decreases  away  from  the  notch,  we 


34 


use  the  finite  element  model  stress  (  oy)  corresponding  to  a 
given  crack  size  in  place  of  the  value  at  the  notch  (  a)  in 
equation  (2.2).  The  edge  crack  solution  is  shown  as  the 
bottom  curve  in  Figure  3.10  with  the  notch  stress  profile 
(oy)  shown  diagrammatically  in  the  inset  figure.  We  expect 
to  see  the  actual  elastic  stress  intensity  solution  for  the 
round  notch  to  approximate  the  edge  crack  solution  for  short 
cracks  and  approach  the  compact  tension  solution  as  the 
crack  size  increases. 

For  the  finite  element  solution,  an  elastic  analysis  was 
performed  for  crack  lengths  up  to  .200  inches  using  the 
finite  element  model.  The  J-integral  was  used  to  evaluate 
the  stress  intensity  solution.  The  J-integral  was  calcu¬ 
lated  over  four  paths  as  shown  in  Figure  3.11.  The 
J  integral  equation  is  given  by  equation  (2.29).  In  the 
elastic  case,  J  is  equal  to  the  energy  release  rate  and  is 
related  to  the  stress  intensity  by  the  following  equation: 

Kj  *  (EJ)l/2  for  plane  stress  (3.5) 

The  stress  intensity  values  calculated  for  the  round  notch 
is  shown  in  Figure  3.12.  These  values  were  averaged  over 
the  J  integral  paths.  Figure  3.12  shows  the  round  notch 
elastic  stress  intensity  to  approximate  the  edge  crack  solu¬ 
tion  for  short  cracks  and  approach  the  compact  tension  spe¬ 
cimen  solution  for  long  cracks.  We  can  also  see  that  the 
notch  influence  on  the  stress  intensity  lies  within  one 
radius  of  the  notch  v approximately  0.2  inches).  This  pheno¬ 
menon  was  also  seen  in  Mercer's  blunt  notch  study.  Thus, 
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given  a  certain  notch  geometry,  we  can  state  that  the  notch 
influence  for  that  particular  geometry  should  lie  within  one 
radius  of  the  notch. 

Crack  Growth  Analysis 

A  crack  growing  from  the  round  notch  of  the  compact  ten¬ 
sion  specimen  was  modeled  using  the  finite  element  model. 
Cracks  up  to  0.024  inches  were  analyzed  under  cyclic  loading 
conditions.  Two  cases  of  loading  were  studied.  The  first 
case  consisted  of  two  load  levels,  P  =  1.148  and  P  =  1.55 
Kips,  and  a  load  ratio  of  0.1.  The  second  case  was  sub¬ 
jected  to  a  plastic  load  level  (P  =  1.55  Kips)  and  a  load 
ratio  of  -1.0.  These  R-ratios  give  a  good  representation  of 
the  stress  levels  which  may  be  seen  in  the  life  cycle  of  a 
component.  These  load  spectra  are  shown  in  Figure  3.13. 

Both  cases  had  a  loading  frequency  of  one  Hz  and  a  constant 
crack  growth  rate.  The  following  results  apply  for  the 
R  »  0.1  case  only. 

Crack  Opening  Displacements 

Figures  3.14  and  3.15  show  the  crack  opening  profiles 
for  crack  lengths  of  0.010,  0.016  and  0.024  inches  at  loads 
of  1.148  Kips  and  1.55  Kips,  respectively.  We  see  that  the 
displacemnts  behind  the  crack  tip  for  the  plastic  load  case 
(Figure  3.15)  are  approximately  1.33  times  greater  than 
those  for  the  elastic  load  case  (Figure  3.14)  for  the  three 
crack  lengths.  Figures  3.16  and  3.17  show  the  crack  pro¬ 
files  at  partial  loads  for  a  crack  which  has  grown  from 
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0.0  inches  to  a  final  length  of  0.024  inches  at  loads  of 
1.148  Kips  and  1.55  Kips,  respectively.  The  cracks  ini¬ 
tially  close  at  the  crack  tip,  but  as  the  load  is  reduced, 
the  crack  closes  at  a  distance  other  than  at  the  crack  tip 
location.  As  the  crack  grows,  a  zone  of  plastic  deformation 
is  left  in  the  wake  of  the  growing  crack.  This  plastic 
deformation  causes  the  shape  of  the  profiles  seen  in  Figures 
3.14-3.17  to  become  concave  downward  and  also  cause  closure 
as  the  load  is  reduced  as  seen  in  Figures  3.16  and  3.17.  We 
again  see  that  even  at  partial  loading,  the  displacements 
behind  the  crack  tip  are  greater  for  the  plar  t '  load  case 
(Figure  3.17)  than  those  seen  for  the  elastic  load  case 
(Figure  3.16).  The  shapes  and  trends  of  Figures  3.14 
through  3.17  are  also  seen  in  Mercer's  (4)  blunt  notch 
study. 

Crack  opening  displacements  at  a  point  .002  inches 
behind  the  crack  tip  were  monitored  as  the  crack  was  grown 
and  are  shown  in  Figure  3.18  for  the  two  different  load 
levels.  The  displacements  shown  are  at  100%  load.  The 
displacements  gradually  increase  as  the  crack  grows  and 
seem  to  be  reaching  a  constant  value  for  the  lower  load 
case.  This  observation  for  the  round  notch  differs  from 
that  seen  by  Mercer  (4).  For  the  blunt  notch,  the  displace¬ 
ments  reached  an  initial  peak  value  early  on  for  the  two 
load  levels,  then  approached  a  constant  displacement  as  the 
crack  continued  to  grow.  Thus,  the  round  notch  geometry 
seems  to  have  an  influence  on  the  results  as  expected. 
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Crack  Closure 


Crack  closure  due  to  the  formation  of  the  plastic  wake 
during  crack  growth  is  seen  in  Figures  3.19  and  3.20  for  the 
two  respective  loads.  ?open  the  opening  load  which 
corresponds  to  the  load  at  which  the  last  node  separates  or 
opens.  Pclose  the  closure  load  which  corresponds  to  the 
load  at  which  the  node  closes  and  exhibits  a  negative 
displacement.  The  closure  and  opening  loads  are  normalized 
with  respect  ;o  the  maximum  load  (PMAX) .  At  the  lower  load 
(Figure  3.19),  we  see  that  the  crack  had  to  grow  at  least 
0.010  inches  before  closure  was  seen  at  that  location.  For 
the  higher  load  (Figure  3.20),  the  crack  had  to  grow  0.008 
inches  before  closure  was  seen  at  that  location.  This 
"delayed"  crack  closure  can  be  attributed  to  the  round  notch 
geometry  condition.  Since  the  round  notch  has  a  lower  notch 
stress  value  than  that  of  a  blunt  notch,  the  associated 
plastic  deformation  develops  slower  than  that  for  the  blunt 
notch.  Thus,  the  crack  must  grow  long  enough,  as  observed 
in  Figures  3.19  and  3.20,  so  that  the  formation  of  a  plastic 
wake  will  develop  and  allow  crack  closure.  Mercer  observed 
some  "delayed"  closure  for  his  elastic  load  case  but  not  for 
his  plastic  load  case.  The  opening  and  closing  loads  gra¬ 
dually  increase  as  the  crack  grows.  This  trend  was  also 
seen  by  Mercer. 

Opening  and  closing  loads  can  be  converted  to  opening 
and  closing  stress  intensities  by  the  following  equations: 
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(3.6) 


Kclose  “  Kmax  ^pclose/pinax^ 

Kopen  =  Kmax  ^popen/pmax^ 

The  values  for  K^x  were  obtained  from  the  previously 
described  elastic  stress  intensity  solution.  Figures  3.21 
and  3.22  show  tho  opening  and  closing  stress  intensities  for 
the  two  load  levels  plotted  versus  the  maximum  stress  inten¬ 
sity.  The  stress  intensities  linearly  increase  with  K^x* 
This  observation  was  also  seen  in  the  blunt  notch  study  (4) 
and  the  single-edge-craeked  work  (5).  One  interesting  point 
regarding  Figures  3.21  and  3.22  is  that  at  the  higher 
applied  load  the  stress  intensity  at  closure  is  lower  than 
that  for  the  elastic  load  at  the  same  Kmax  level.  This  phe¬ 
nomenon  was  also  seen  in  the  single -edge-cracked  specimen, 
as  shown  in  Figure  3.23,  and  the  blunt  notch  crack  closure 
studies.  Thus,  if  a  parameter  such  as  fiKeff  *  Kmax  “  Kclose 
is  a  driver  of  crack  growth,  the  higher  load  would  result  in 
a  larger  effective  stress  intensity  range  than  the  lower 
load  at  the  same  maximum  applied  stress  intensity. 

Stress  Intensity  Ranges 

As  stated  earlier,  the  crack  growth  rate  may  be  governed 
by  a  4Kejs£  parameter.  Thus,  we  need  to  evaluate  the 
influence  of  the  Stress  level  on  the  growth  rate  of  the 
round  notch  crack.  This  is  best  done  by  evaluating  the 
stress  intensity  ranges  for  the  two  applied  loads. 

For  an  applied  load  of  1.148  Kips,  Figure  3.24  shows  the 
maximum  stress  intensity,  minimum  stress  intensity,  and 
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stress  intensity  at  closure  for  a  crack  grown  to  0.024 
inches.  The  minimum  stress  intensity  is  taken  as  10%  of  the 
maximum  stress  intensity.  From  this  figure  we  can  define  a 
corresponding  stress  intensity  range  (  A K)  and  an  effective 
stress  intensity  range  (  &Keff}  by  tfie  following  equations: 

AK  =  Kffiax  -  Km£n  (3.8) 

iKeff  *  Kraax  ”  ^close  (°r  Kopen*  (3.9) 

These  ranges  are  shown  in  Figure  3? 25.  We  see  that  over 

the  0.024-inch  crack  growth  range,  AK  increases  by  124%, 
but  AKeff  changes  only  19%. 

At  the  higher  load  of  1.55  Kips,  we  see  the  values  of 
Kmax  '  Kmin'  and  Kclose  *n  Figure  3.26.  Figure  3.27  shows 
the  corresponding  AK  and  AKeff  values.  In  Figure  3.27, 

K  increases  by  123%,  while  AKeff  changes  by  22%. 

We  observe  that  AKeff  is  greater  at  the  higher  load 
level  (Figure  3.27)  since  Kciose  (or  Kopen)  is  smaller  at  a 
common  value  of  K^x  of  approximately  32  Ks£  /in  as  seen  in 
Figures  3.21  and  3.22.  Thus,  we  can  say  that  the  crack 
effect  at  the  round  notch  should  be  more  pronounced  at  the 
higher  load  level  than  the  lower  one.  This  observation  was 
also  seen  in  the  aingle-edge-crack  work  (5)  and  can  be  seen 
in  Figure  3.23. 

Stresses  and  Strains 

The  stresses  and  strains  at  and  ahead  of  the  crack  tip 
were  monitored.  The  values  at  the  crack  tip  actually  refer 
to  the  values  calculated  at  the  element  integration  point 


closest  to  the  crack  tip.  For  the  eight-noded  quadratic 
elements,  this  point  is  located  approximately  0.0004  inches 
from  the  crack  tip. 


The  stress  response  at  the  crack  tip  ranged  between  full 
tensile  and  compressive  yielding  for  both  load  cases  and  is 
shown  in  Figure  3.28  for  the  higher  load  case.  This  results 
in  a  crack  tip  R  ratio  of  approximately  -1.0. 

The  crack  tip  strains  were  found  to  increase  similarly 
to  that  seen  for  the  maximum  elastic  stress  intensity  as 
shown  in  Figures  3.29  and  3.30  fovr  t.he  two  load  cases.  As 
expected,  the  maximum  strains  are  g.vettnr  for  the  higher 
load  level  than  for  the  lower  load  Level.  These  trends  were 
also  seen  by  Mercer  (4). 

The  stress  profiles  ahead  of  the  crack  tip  for  crack 
lengths  of  .010,  .016,  and  .024  inches  are  shown  in 
Figures  3.31-3.34  at  100%  and  10%  loading  for  the  respective 
load  cases.  In  Figure  3.31,  we  see  that  the  stress  levels 
remain  constant  for  a  short  period  near  the  notch,  but 
steadily  decrease  as  the  crack  grows  away  from  the  notch. 

For  the  10%  load  condition  (Figure  3.32) ,  we  see  that  the 
stresses  dramatically  increase  to  a  peak  value  and  then 
gradually  decrease  to  a  constant  value.  Tt  •  plastic  zone 
extent  is  shown  in  Figure  3.32  and  is  approximately  the  same 
size  as  the  three  respective  crack  sizes.  This  is  typical 
of  short  crack  behavior.  At  the  higher  load  level  (Figure 
3,33),  we  see  that  the  stress  levels  remain  constant  similar 
to  Figure  3.31  but  for  a  longer  duration.  The  stresses 
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remain  constant  as  the  crack  grows  within  the  notch  plastic 
zone  and  starts  to  decrease  once  the  crack  grows  beyond  this 
plastic  region.  In  Figure  3.34,  we  again  see  that  the  crack 
tip  plastic  zone  extent  for  the  three  crack  lengths  are 
approximately  the  same  size  as  the  cracks.  Figure  3.34  is 
similar  in  shape  to  Figure  3.32,  but  the  plastic  zone  extent 
is  more  pronounced  in  Figure  3.34  due  to  the  plastic  load 
condition. 

The  strain  profiles  ahead  of  the  crack  tip  for  crack 
lengths  of  .010,  .016,  and  .024  inches  are  shown  in 
Figures  3.35-3.38  at  100%  and  10%  loading  for  the  two  load 
cases.  At  the  lower  load  (Figure  3.35),  we  see  that  the 
strains  are  greatest  at  the  crack  tip  and  gradually 
decrease  away  from  the  notch.  As  the  strain  decreases,  it 
approaches  a  constant  value  for  the  three  different  crack 
sizes.  The  plastic  zone  extent  is  shown  in  Figure  3.35  and 
is  approximate] y  the  same  size  as  the  respective  crack 
length.  As  stated  earlier,  this  phenomenon  is  typical  of 
short  crack  behavior.  Figure  3.36  shows  the  strain  profile 
ahead  of  the  crack  tip  at  10%  load  for  the  elastic  load 
case.  The  shape  is  similar  to  Figure  3.35  but  at  lower 
strain  levels.  At  the  higher  load  level  (Figure  3.37),  we 
see  the  expected  result  of  the  strain  values  being  greater 
for  the  higher  load  than  for  the  lower  load.  As  the  crack 
grows  out  of  the  crack  tip  plastic  2one,  we  see  that  the 
strain  values  are  approximately  the  same  for  the  high  load 
case  and  the  low  load  case,  We  also  see  the  same  phenomenon 
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noticed  in  Figure  3.35  regarding  the  plastic  zone  extent. 
Figure  3.38  shows  the  strain  profile  at  10%  load  for  the 
plastic  load  case.  The  curve  is  again  similar  in  shape  to 
the  previous  strain  profiles.  The  strain  values  are 
slightly  higher  than  those  seen  in  Figure  3.36,  as  expected. 

A  comparison  of  plastic  strain  profiles  can  be  seen  from 
Figures  3.39  and  3.40  for  the  two  load  cases.  The  higher 
load  induces  significantly  more  plasticity  than  the  lower 
load.  For  a  crack  length  of  .024  inches,  Figure  3.41  shows 
the  accumulated  plastic  strain  ahead  of,  as  well  as  behind, 
the  crack  tip  is  greater  for  the  higher  load  case  than  the 
lower  load  case.  As  the  crack  tip  grows,  a  region  of 
plasticity  or  plastic  wake  is  formed.  This  region  is  pro¬ 
duced  between  the  original  crack  tip  location  and  its 
current  position.  Ahead  of  the  crack  tip,  we  see  from 
Figure  3.41  that  the  amount  of  plastic  strain  appears  to  be 
less  than  that  behind  the  crack  tip.  Thus,  the  formation  of 
the  plastic  wake  along  the  crack  surface  seems  to  restrict 
the  amount  of  plasticity  ahead  of  the  crack  tip.  In  Figure 
3.41,  the  plastic  wake  increases  dramatically  behind  the 
crack  tip  for  the  two  load  level  cases  and  the  amount  of 
plastic  strain  seems  to  vary  directly  in  relation  to  the 
magnitude  of  the  load  applied.  The  plastic  zone  size  for 
the  .024-inch  crack  is  seen  to  be  approximately  .025  inches 

for  plastic  strains  greater  than  0.5%.  The  spikes  observed 
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in  the  curves  are  attributed  to  numerical  inaccuracies 
within  the  finite  element  program  due  to  a  plastic  strain 


singularity  effect  at  the  crack  tip;  at  this  point  the 
theory  of  small  displacement  is  no  longer  valid. 

The  stress  and  strain  trends  observed  for  cracks  growing 
from  a  round  notch  specimen  were  similar  to  those  seen  in 
the  blunt  notch  work  (as  shown  in  Figures  3.42  and  3.43)  and 
the  single-edge-crack  specimen  (as  shown  in  Figure  3.44). 
Both  of  these  studies  centered  on  short  crack  growth  analy¬ 
sis. 

Negative  R-Ratio  Results 

The  results  and  discussion  in  this  section  are  derived 
from  the  compact  tension  specimen  loaded  at  the  higher  load 
level  and  with  a  load  ratio  of  -1.0.  The  results  for  the 
negative  load  ratio  can  be  broken  down  into  three  major 
areas : 

a.  Crack  opening  displacements  at  and  behind  the  crack 

tip. 

b.  Opening  and  closure  load  information  during  crack 
growth . 

c.  Profiles  of  stress  and  plastic  strain  ahead  of  the 
crack  tip. 

Crack  opening  profiles  for  crack  lengths  of  .010,  .016, 
and  .024  inches  are  shown  in  Figure  3.45  for  the  higher  load 
case.  The  curves  for  the  negative  load  ratio  case  are  simi¬ 
lar  to  those  for  the  positive  load  ratio  case  (Figure  3.15), 
but  we  see  that  the  displacements  behind  the  crack  tip  for 
the  three  crack  lengths  are  approximately  1.5  times  greater 
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for  the  negative  load  ratio  than  for  the  positive  load  ratio 
condition. 

Crack  opening  displacements  at  a  point  .002  inches 
behind  the  crack  tip  are  observed  in  Figure  3.46.  These 
displacements  are  approximately  three  times  greater  than 
those  seen  in  the  positive  load  ratio  case  (Figure  3.18). 

Crack  opening  and  closure  loads  were  examined  for  a 
point  .002  inches  behind  the  crack  tip.  These  loads  are 
shown  in  Figure  3.47  for  which  crack  closure  and  opening 
occurred  for  the  negative  load  ratio  condition.  We  see  that 
for  the  plastic  load  the  crack  closes  immediately.  However, 
for  the  positive  load  ratio  case  with  the  same  loading 
(Figure  3.20),  the  crack  had  to  grow  .008  inches  before  clo¬ 
sure  was  seen  at  that  location.  Thus,  crack  closure  devel¬ 
ops  quicker  for  R  *  -1  than  for  R  «  0.1  for  the  same  load 
value.  The  closure  loads  decrease  gradually  with  increasing 
crack  extension  and  the  opening  loads  approach  a  constant 
value  as  the  crack  grows. 

Opening  and  closure  stress  intensities  are  shown  in 
Figure  3.48  for  the  plastic  loading  case.  The  closing  and 
opening  stress  intensities  are  defined  by  equations  (3.6) 
and  (3.7),  respectively.  The  values  for  the  opening  stress 
intensity  remain  relatively  constant  while  the  closure 
stress  intensity  values  decrease  steadily  with  increasing 
crack  growth.  These  opening  and  closure  load  trends,  for  a 
crack  growing  from  a  round  notch  with  a  negative  load  ratio, 
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are  similar  to  those  seen  in  the  single -edge -cracked  speci¬ 
men  work  (5),  as  shown  in  Figure  3.49. 

The  stress  profiles  ahead  of  the  crack  tip  at  100%  and 
-100%  load  for  crack  lengths  of  .010,  .016,  and  .024  inches 
are  shown  in  Figures  3.50  and  3.51,  respectively.  For  the 
negative  load  ratio  case  at  100%  load  (Figure  3.50),  we  see 
results  that  are  similar  to  those  found  in  the  positive  load 
ratio  case  at  100%  load  (Figure  3.33).  At  the  full 
compressive  load  (Figure  3.51),  we  see  the  stress  values 
remain  relatively  constant,  then  gradually  decrease. 

The  plastic  strain  profile  ahead  of  the  crack  tip  for 
the  negative  load  ratio  case  is  seen  in  Figure  3.52.  The 
accumulated  plastic  strain  ahead  of,  as  well  as  behind,  the 
crack  tip  for  a  crack  length  of  .024  inches  is  shown  in 
Figure  3.53  for  the  plastic  load  case.  We  see  that  the 
plastic  strains  behind  and  ahead  of  the  crack  tip  are 
smaller  for  the  negative  load  ratio  than  those  seen  for  the 
positive  load  ratio  case  (Figure  3.41)  at  the  same  load 
(P  *=  1.55  Kips).  This  indicates  that  there  is  a  severe 
compression  developed  in  the  plastic  wake  under  fully 
reversed  loading  and  these  compressive  forces  have 
restricted  the  size  of  the  plastic  wake  in  the  region 
directly  behind  the  crack  tip.  This  phenomenon  was  also 
observed  in  the  single-edge-cracked  work  (5),  as  seen  in 
Figure  3.44.  The  spike  seen  at  the  crack  tip  is  attributed 
to  numerical  inaccuracies  within  the  program,  as  described 
earlier. 
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Figure  3.2:  Cracked  Notched  Specimen  Overall  Model 
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Plastic  Zone  Size  versus  Applied  Load,  Viscoplastic  Analysis 
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Uncracked  Round  Notch  Stress-Strain  Response  During  Cyclic 
Loading 


Round  Notch  Stress  Intensity  Solution  Bounds 


Figure  3.11s  J-Integral  Paths 


Figure  3.12:  Stress  Intensity  Solution  for  Round  Notch  using  J-Integral 


Figure  3.13t  Load  Ratio  Cases  Used  in  Round  Notch  Study 
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Figure  3.14:  Crack  Opening  Profile,  p=  1.148  Kips,  R=  0.1,  10096  Load 
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Figure  3.15:  Crack  Opening  Profile,  P=  1.55  Kips,  R=  0.1,  100%  Load 


Figure  3.16:  Crack  Opening  Profile  at  Partial  Loading,  a=  .024",  P=  1.148  Kips 
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Figure  3.17:  Crack  Opening  Profile  at  Partial  Loading,  a=  .024",  P=  1.55  Kips, 
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Figure  3.18:  Crack  Opening  Displacement  versus  Crack  Size  for  P=  1.148  and  P=  1.55 


Figure  3.21:  Closure  .002"  Behind  Crack  Tip  versus  Maximum  Stress  Intensity, 
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Figure  3.22:  Closure  .002"  Behind  Crack  Tip  versus  Maximum  Stress  Intensity, 


Specimen  (Ref.  5) 


Figure  3.24:  Stress  Intensity  versus  Crack  Size,  P=  1.148  Kips 
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Figure  3.25:  Stress  Intensity  Ranges  versus  Crack  Size,  P=  1.148  Kips 


Figure  3.26:  Stress  Intensity  versus  Crack  Size,  P=  1.55  Kips 


Figure  3.27:  Stress  Intensity  Ranges  versus  Crack  Size,  P=  1.55  Kips, 
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Figure  3.28:  Crack  Tip  Stresses,  P=  1.55  Kips 
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Figure  3.31:  Stress  Profile  Ahead  of  Crack,  P=  1.148  Kips,  R=  0.1,  1009S  Load 
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Figure  3.32:  Stress  Profile  Ahead  of  Crack,  P=  1.148  Kips,  R=  0.1,  10%  Load 


Figure  3.33:  Stress  Profile  Ahead  of  Crack,  P=  1.55  Kips,  R=  0.1,  100%  Load 
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.148  Kips,  R=  0.1,  10%  Load 
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Figure  3.37:  Strain  Profile  Ahead  of  Crack,  P=  1.55  Kips,  R=  0.1,  100%  Load 
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Figure  3.38:  Strain  Profile  Ahead  of  Crack,  P=  1.55  Kips,  R=  0.1,  10%  Load 


Figure  3.39:  Plastic  Strain  Profile  Ahead  of  Crack,  P=  1.55  Kips 


10*01 


* 


♦  • 
O  *4 


II  II 
06  CL 


86 


Figure  3.40:  Plastic  Strain  Profile  Ahead  of  Crack,  P=  1.148  Kips,  R=  0.1,  100%  Load 
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Figure  3.41:  Accumulated  Plastic  Strain,  a=  .024",  P=  1.148  and  P=  1.55  Kips, 
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Figure  3.43:  Accumulated  Plastic  strain  for  Blunt 
Notch  Specimen  (Ref.  4) 
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Figure  3.45:  Crack  Opening  Profile,  P=  1.55  Kips,  R=  -1.0,  100%  Load 
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Figure  3.47:  Closure  and  Opening 
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Figure  3.48:  Closure  and  Opening  Stress  Intensities  versus  Maximum  Stress  Intensity 
for  P=  1.55  Kips,  R=  -1.0 
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Figure  3-49:  Crack  Closure  Stress  Intensities  for  Single-Edge-Cracked  Specimen 
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Figure  3.52:  Plastic  Strain  Profile  Ahead  of  Crack  Tip,  P=  1.55  Kips 
100%  Load 
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Figure  3.53:  Accumulated  Plastic  Strain,  a=  .024",  P=  1.55  Kips,  R= 


IV.  Summary  and  Conclusions 


A  two-dimensional  finite  element  code  called  SNAP  was 
used  to  model  crack  growth  from  a  round  notch  compact  ten¬ 
sion  specimen  of  IN  718  at  1200°.  The  nonlinear  material 
behavior  was  modeled  using  the  Bodner-Partom  viscoplastic 
flow  law  and  an  implicit  scheme  was  used  to  integrate  these 
plastic  strain  relations  with  respect  to  time. 

The  IN  718  specimen  was  modeled  under  cyclic  conditions 
at  two  load  levels  with  a  load  ratio  of  0.1.  The  specimen 
was  also  analyzed  under  plastic  loading  with  a  load  ratio  of 
-1.0. 

The  following  statements  and  conclusions  are  based  on 
the  fatigue  analysis  presented  herein: 

1.  For  cyclic  loading  frequencies  of  1.0  Hz,  it  was 
found  that  the  majority  of  plastic  strain  at  the  crack  tip 
is  accumulated  in  the  first  load  cycle. 

2.  For  the  round  notch  specimen,  an  elastic  stress  con¬ 
centration  factor  of  2.18  was  obtained  numerically.  This 
result  compared  favorably  to  that  found  by  Wilson  (24)  and 
was  approximately  one-half  the  value  found  by  Mercer  (4)  for 
a  blunt  notch  specimen.  Thus,  the  round  notch  geometry  was 
seen  to  have  lower  notch  stress  values  than  that  of  a  blunt 
notch  geometry. 

3.  An  elastic  stress  intensity  solution  for  a  crack 
growing  from  the  round  notch  was  determined  using  the  SNAP 
finite  element  model.  The  round  notch  stress  intensity  was 
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seen  to  approximate  the  edge  crack  solution  for  short  cracks 
and  approach  the  compact  tension  specimen  solution  for  long 
cracks.  The  notch  influence  region  on  the  elastic  stress 
intensity  was  observed  to  be  on  the  order  of  one  notch 
radius  (approximately  0.2  inches).  This  phenomenon  was  also 
seen  by  Mercer  in  his  blunt  notch  work. 

4.  In  general,  the  crack  growth  trends  (i.e.,  stress- 
strain  values,  crack  opening  displacements,  and  opening  and 
closure  information)  for  a  crack  growing  from  a  round  notch 
agree  with  those  seen  in  the  blunt  notch  and  single-edge- 
cracked  specimens.  Thus,  we  can  optimistically  state  that 
these  trends  observed  are  those  typically  encountered  in 
short  crack  studies.  This  short  crack  behavior  is  easily 
seen  from  the  stress  and  strain  values  ahead  of  the  crack 
tip,  where  the  crack  tip  plastic  tone  is  approximately  the 
same  size  as  the  crack  length. 

5.  *'or  the  load  ratio  case  of  Q.l,  the  crack  had  to 
grow  a  significant  amount  before  closure  was  observed.  This 
"delayed"  crack  closure  is  attributed  to  the  round  notch 
geometry  condition.  Por  the  load  ratio  case  of  -1.0,  crack 
closure  was  immediately  encountered.  Thus,  closure  develops 
quicker  for  a  negative  load  rrvio  condition  than  for  a  posi¬ 
tive  load  ratio  case.  This  observation  was  also  seen  in  the 
single-edge-cracked  specimen  (5). 

It  is  believed  that  the  round  notch  fatigue  crack 
results  presented  herein  will  add  to  the  body  of  knowledge 
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concerning  high  temperature  cyclic  crack  growth  and  may  pro¬ 


vide  some 
growth  of 


insight  into  the  problems  associated  with  fatigue 
short  cracks. 
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Crack  closure  was  seen  to  develop  overa~Tbnger  distance 
for  the  round  notch  than  that  for  a  blunt  notch.  The  elastic 
stress  intensity  influence  on  the  round  notch  was  found  to 
be  within  one  radius  of  the  notch.  In  general,  the  crack 
growth  trends  observed  for  the  round  notch  are  similar  to 
those  seen  in  the  blunt  notch  and  single-edge-cracked  studies. 


